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Abstract 

This report documents the procedures developed for incorporating smart laminate and panel analysis 
capabilities within the HyperSizer aerospace structural sizing software package. HyperSizer analyzes 
stiffened panels composed of arbitrary composite laminates through stiffener homogenization, or 
“smearing”, techniques. The result is an effective constitutive equation for the stiffened panel that is 
suitable for use in a full vehicle-scale finite element analysis via MSC/NASTRAN. The existing thermo- 
elastic capabilities of HyperSizer have herein been extended to include coupled thermo-electro-magneto- 
elastic analysis capabilities. This represents a significant step toward realization of design tools capable of 
guiding the development of the next generation of smart aerospace structures. Verification results are 
presented that compare the developed smart HyperSizer capability with an ABAQUS piezoelectric finite 
element solution for a facesheet-flange combination. These results show good agreement between 
HyperSizer and ABAQUS, but highlight a limitation of the HyperSizer formulation in that constant 
electric field components are assumed. 


1. Introduction 

Adaptive structures show a great deal of promise for future aerospace applications. The envisioned 
structure’s adaptive capabilities will rely on so-called “smart” (or “intelligent”) materials, which have 
properties that enable them to sense various stimuli and react in some way. Incorporating these intelligent 
materials has the potential to remove fundamental design constraints and transform aerospace structures 
into life-like responsive systems, enabling optimum performance, reliability, and weight throughout a 
changing mission profile. The work presented herein focuses on the structural application of one type of 
smart material: piezo-electro-magnetic materials (i.e., piezoelectric and piezomagnetic ceramics). Piezo- 
electro-magnetic materials are those that exhibit coupling among their electric, magnetic, mechanical, and 
thermal responses. That is, for example, in response to an applied voltage or current, a piezoelectric 
material (such as lead zirconium titanate, or PZT) will respond mechanically with a change in strain or 
stress. Conversely, if piezoelectric materials are loaded mechanically or thermally, a change in their 
electric field or flux results. Thus, piezoelectric materials can serve as both actuators and sensors and have 
the potential to perform tasks passively through utilization of an induced field/flux to cause a desired 
mechanical response. Piezomagnetic materials (such as CoFe 2 0 4 ) exhibit similar coupling between then- 
magnetic and mechanical behaviors, and by forming a composite of piezoelectric and piezomagnetic 
materials, a fully coupled piezo-electro-magnetic material can be produced. Piezo-electro-magnetic 
materials are characterized by their fast response times to applied (or sensed) stimuli (on the order of 10“ 2 
to 10 1 milliseconds) and thus have found significant application as vibration dampers. Reviews of piezo- 
electro-magnetic concepts and materials are available in references 1 to 4. 

In order to realize the potential embodied by smart materials and structures, advances in modeling 
and simulation technologies are needed. The standard tools for structural design are finite element 
analysis (FEA) models (e.g., ANSYS, NASTRAN, ABAQUS). However, FEA models are ill-suited 
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(i.e., inefficient, subject to operator error) for rapid design and sizing (i.e., trade studies) for structural 
components. Further, the lack of well-developed and robust capabilities related to intelligent materials 
underscores the shortcomings of the FEA approach when it comes to adaptive structures. There is thus a 
need for physics-based design, analysis, and sizing tools are needed that capture the essential 
characteristics of piezo-electro-magnetic materials and enable the analysis of structures composed of 
these materials. 

Major investments have been made by NASA over the last decade that laid the groundwork for such 
design tools. The work described herein has resulted from the research project entitled “Multi-Scale 
Sizing of Lightweight Multifunctional Spacecraft Structural Components.” This project, funded by NASA 
Headquarters, has brought together NASA Glenn Research Center’s Micromechanical Analysis Code 
with Generalized Method of Cells (MAC/GMC) (refs. 5 and 6), which simulates the nonlinear behavior of 
smart and composite materials, and HyperSizer (ref. 7), a commercial structural sizing software package 
originating from NASA Langley technology. Both MAC/GMC and HyperSizer have been enhanced to 
simulate piezo-electro-magnetic materials and seamlessly linked such that MAC/GMC provides the ply- 
level behavior of a traditional or smart composite (or monolithic) within a stiffened structure modeled by 
HyperSizer. In addition, HyperSizer has been enhanced to enable consideration of time-dependent 
loading, allowing simulation of an entire mission profile. As such, the software can now consider many 
points from a structure’s operating envelope, rather than simply sizing based on a single load level. 
HyperSizer also links with FEA to enable automatic application of higher-scale structural loads on the 
structural components that are optimized by the software. Thus the integrated HyperSizer — MAC/GMC 
product now represents a unique multi-scale tool for the analysis of advance lightweight aerospace 
structures. Through its linkage with FEA, the software can consider a truly integrated vehicle structural 
design rather than an isolated design of each component. 

This report describes the methods and procedures that have been developed to enable the analysis of 
piezo-electro-magnetic materials within HyperSizer. Starting with the analysis of a piezo-electro- 
magnetic laminate, new laminate level matrices that account for the electric, magnetic, thermo-electric 
and thermo-magnetic effects are developed. These then can be treated in a way analogous to the 
HyperSizer treatment of the laminate thermal matrix, enabling use of homogenization, or “smearing” 
techniques to develop stiffened panel level electric, magnetic, thermo-electric, and thermo-magnetic 
terms. The developed methods are based on a classical lamination theory treatment of the laminates 
comprising a given stiffened panel, and the homogenization of the stiffened panel so it can be represented 
with classical lamination theory terms. It is this simplicity that provides the methods with the level of 
efficiency needed to consider many design cases rapidly while still capturing the dominant first-order 
effects. A good body of work exists for the analysis of piezoelectric laminates, using both analytical 
(refs. 8 to 17) and finite element (refs. 18 to 22) approaches. The piezomagnetic laminates have also 
received some attention (refs. 23 to 25). The work most closely related to the methods developed herein 
involve the extension of classical lamination theory to include piezoelectric plies by Lee (ref. 8), Crawley 
and Lazarus (ref. 9) and Tauchert (ref. 10), and to include piezomagnetic and inelastic behavior by 
Bednarcyk (ref. 26). 


2. Reference Plane Shifting Procedure for a 
Thermo-Electro-Magneto-Elastic Laminate 

Because a stiffened panel is, in general, composed of a number of laminates, each of which has 
properties that are typically calculated with respect to its own midplane, a first step in calculating 
properties of the panel as a whole involves shifting the reference plane of a given laminate. As will be 
shown, if the reference plane of each laminate comprising a stiffened panel can be shifted to a common 
reference plane, the homogenization of the laminate properties to form the properties of the stiffened 
panel properties becomes straightforward. 
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Consider and arbitrary number of composite laminates, each of which has its constitutive behavior 
defined with respect to its midplane via the standard thermo-elastic lamination theory equation (refs. 27 
and 28), 
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In this equation, N and M are the force and moment resultant vectors, which are related to the midplane 
strain and curvature vectors, e° and K°, by the laminate extensional, coupling, and bending stiffness 
matrices (measured with respect to the midplane), A 0 , B°, and D°, and the thermal force and moment 

resultant vectors (measured with respect to the midplane), Nq and Mq . The laminate midplane stiffness 
matrices and thermal resultants are given by, 
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where is the reduced stiffness matrix of ply k, a* is the coefficient of thermal expansion vector of ply 

k, z\ is the z-coordinate position of the top of ply k measured with respect to the laminate midplane, A To 
is the temperature change from reference temperature at the laminate midplane, and AG is the linear 
laminate through-thickness temperature gradient, i.e., T(z) - T Kf = A T 0 +z AG, with T(z) being the 
temperature at any z-coordinate location and being the reference temperature. 

We now consider a shift of the reference plane from the laminate midplane to an arbitrary z- 
coordinate position and seek to determine the effect on the quantities given in eqs. (2) and (3). This 
arbitrary reference plane shift is depicted in figure 1 . Considering point A, in the original, laminate 

midplane coordinate system, this point is located at z® =-/z/2 , whereas, in the new shifted coordinate 
system, this point is located at z^ ew = - h/2 - Az . Thus, for an arbitrary shift of reference plane we have 
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where Az measures the distance in the positive z-direction from original laminate reference plane to the 
new laminate reference plane. The temperature change from reference temperature at the new reference 
plane is given based on the linear through-thickness temperature gradient as, 

AT new =AT 0 +AGAz (5) 

The laminate stiffness matrices, measured with respect to the new shifted reference plane are given 

by, 
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Figure 1 . — Arbitrary reference plane shift from the laminate midplane. 
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Substituting for z new in eqs. (6) to (8) using eq. (4) we arrive at, 
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Likewise, for the thermal force and moment resultants measured with respect to the new shifted 
reference plane, we have, 
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Substituting for z new in eq. (12) using eq. (4) and for A T mw using eq. (5), we arrive at, 
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HyperSizer employs an alternate non-classical form of the laminate constitutive eq. (1) in which the 
thermal effects are accounted for using “thermal ABD” terms (ref. 29). In this case, the laminate 
constitutive equation is written as, 
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the thermal force and moment resultants are related to the thermal ABD terms by, 
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With respect to the new shifted reference plane, we have, 
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Substituting eq. (5) into eq. (18) yields, 
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Using eq. (13) and comparing the terms present in eq. (20) with those in eq. (16) gives, 
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Comparing the terms in eq. (25) with those in eq. (24), we confirm eq. (23) and also have, 


Dnew - Bnew Az = D“ + Bq Az (26) 

Rearranging eq. (26) and substituting using eq. (23) then yields, 

Dnew = Dq + 2AzB“ + Az 2 A“ (27) 


Equations (9) to (1 1), (13), (14), (21), (23), and (27) thus enable the determination of the thermo- 
elastic laminate constitutive equation terms with respect to an arbitrary reference plane from the reference 
plane thermo-elastic laminate constitutive equation terms and the reference plane shift, Az . 

Thermo-piezo-electro-magnetic effects can be included within the midplane laminate constitutive 
eq. (1) (see ref. 26 for details) as, 
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where N 0 and M 0 are the electric force an moment resultant vectors, N 0 and M 0 are the thermo- 
electric force and moment resultant vectors, and are the magnetic force and moment resultant 

1 jrw I 1 jfTTl 

vectors, and Nq and M” are the thermo-magnetic force and moment resultant vectors, all determined 
with respect to the laminate midplane. These thermo-piezo-electro-magnetic force and moment resultant 
vectors are given by, 
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where E t are the electric field components, Ej are the thermo-electric field components, H t are the 

magnetic field components, Hf are the thermo-magnetic field components, e k is the 3x3 reduced 
piezoelectric coefficient matrix (for layer k ), 
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and is the 3x3 reduced piezomagnetic matrix (for layer k). 
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(see ref. 26 for details). It should be noted that the minus signs present in the above moment resultant eqs. 
(29) to (32) are due to the coordinate system employed shown in figure 1. This coordinate system differs 
from that used by Bednarcyk (ref. 26), in which the aforementioned minus signs do not appear. 

For consistency, we now develop an alternative representation of these force and moment resultants 
that by introducing terms analogous to the thermal ABD terms present in eq. (15). Restricting the applied 
electric field components to the case in which they are constant throughout the laminate, eq. (29) can be 
written as, 
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Combining eq. (33), 
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or, 
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where, 


A o = Z & ] ( z k-\ ~ z k) B o = “Z [®* 1 ( z k-\ ~ z l ) 


Here Af and Bf are each 3x3 matrices and combine to form the laminate electric AB matrix, as shown 
in eq. (35). 

Similarly, restricting the applied magnetic field components to the case in which they are constant 
throughout the laminate, we can write eq. (31) as, 
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Here A™ and are each 3x3 matrices and combine to form the laminate magnetic AB matrix, as 
shown in eq. (37). 

The temperature change from reference at any point within a ply is given by T(z ) - T ret = ATq + z AG, 
and the thermo-electric and thermo-magnetic field components, appearing in eqs. (30) and (32), are 
related to this temperature change by, 


El =C k {AT 0 +z AG) 


Hi =M4(Ar 0 +zAG) 






where ^ is the effective pyroelectric constant vector of ply k and \\i* k is the effective pyromagnetic 
constant vector for ply k, both in the global laminate coordinate system. Note that the relations between 
the thermal field quantities and the pyroelectric and pyromagnetic constants in eq. (39) is analogous to 
that of the thermal strain and the coefficient of thermal expansion. Substituting eq. (39) into eqs. (30) and 
(32) yields, 


NASA/CR— 2009-21 5269 


11 



N z k 


N z k 


No 7 =Efrt] | £ (A T 0 +z A G)dz Mf =-'£[e k ] { £ (A7b +z AG)zcfe (40) 


AT z * 


AT z * 


No® 1 =£[«*] I ¥i( A7 0 +*AG)& M^=-£[q*] j y* k (AT 0 +z AG)zdz (41) 


Nfl ; 
Mf j' 


Nfl. 
Mo®" J 


X [®* ] £ ( z k-\ - z k ) AT’o + tZ [** ] £ (4-1 “ 4 ) AG 

*=1 z k=l 

~~Z [®* ] £ (4-1 - 4 ) A7o - “X [®* ] £ (4-1 1 “ 4 ) AG 
. z Jt=l 3 k = 1 J 

Z [fl* ] ¥* ( z Jt-l - z * ) AT o + ^Z [4* 1 V* (4-1 - 4 ) AG 

*=1 1 k = l 

-^ZfaM (4-i -4)Ar 0 -IZ[q*K(4-i -4)ag 


Nfl = rAf BflrATol 
Mfj [ B o r Pg r jbgJ 


N^ 7 1_ A^ 7 B^lr A^o ] 
M^ 7 J B^ 7 D^ 7 J _-Ag[j 


Af =£[«*]£ (*«-**) B f r =“Z[®*l£(4-i-4) D o r =^ZK]4(4-,-4) ( 45 ) 


A o /r= Z [q*]v*( z *-i- z *) B r=4Z[q*M(4-i-4) D <T=^Z[q*]v*(4-i-4) ( 46 ) 


Note that each of the above thermo-electric and thermo-magnetic ABD terms is a 3x1 matrix. We can 
now write the laminate constitutive eqs. (15) and (28) as, 


NASA/CR— 2009-21 5269 


12 



( 47 ) 
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A7 o,l 
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'Ar 0 1 
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-AG] 

j>MT 

L^o 

d"J 


Comparing the electric, magnetic, thermo-electric, and thermo-magnetic terms in eq. (47) with the 
thermal terms in eq. (47) (see eqs. (15a) to (15c), (36), (38), and (45) to (46), it is clear that there is a 

direct analogy among all A J 0 matrices, among all matrices, and among all matrices. Since the z k 

functionality in all of these thermo-electro-magnetic ABD terms is the same as that of the thermal ABD 
terms, it can be shown that these thermo-electro-magnetic ABD terms will shift reference plane (see 
fig. 1) in a way analogous to the shift in the thermal ABD terms given in eqs. (21), (23), and (27). Thus, 
for a reference plane shift of Az , we have, 

A4 W =A^ BJ new =B J 0 +AzA J 0 <, j = E, M, ET, MT (48) 

D 4 W = D' + 2AzB' + Az 2 A^ , j = ET, MT (49) 

In terms of the thermo-electro-magnetic force and moment resultants (eqs. (33), (38), (42), and (43)), we 
have, 


Nnew = Nq , j = E, M, ET, MT 

(50) 

=M^ + AzN' , j = E, M, ET, MT 

(51) 


A simple procedure has thus been established for shifting all elements of the thermo-electro-magneto- 
elastic laminate constitutive eq. (47) from the laminate midplane to an arbitrary reference plane. This new 
laminate constitutive equation is given by, 


Nl_ 

^new jjnew 

£ neW 

1 ' 

M] 

jjnew j^new j 

_K neW 

J . 


new 
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n ET i 
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E y 

E z'\ 


a m 1 

^new 1 
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H y 

V H z'\ 


J) MT 1 1 
"new U 


new 

-AG ] 


(52) 


3.0 Homogenization Procedure for a Group of Laminates 

Now that a reference plane shifting procedure has been established, it is possible to develop a simple 
methodology for homogenizing the constitutive equations of a group of laminates to arrive at an effective 
laminate constitutive equation of the laminate constituted by the group. Because each term in the laminate 
constitutive eq. (47) is determined via summation of ply-level quantities, provided that the constitutive 
equations of all laminates in the group are determined with respect to the same reference plane (eq. (52)), 
the constitutive equation of the group can be determined via a simple summation. Thus, given a number 
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Z NEW 



Figure 2. — A group of laminates whose constitutive equation terms have been 
determined with respect to a common reference plane, defined by the origin of 
the coordinate z ne w- This allows determination of effective constitutive equation 
terms (with respect to the common reference plane) for the laminate constituted 
by the group. 


of laminates and the z new coordinate of the original reference plane of each laminate (given by l h see 
fig. 2), we can shift the reference plane of each laminate to the common z new reference plane and then sum 
the resulting ABDp matrices to determine the new homogenized ABD of the entire group of laminates. 

For each laminate (denoted by the subscript p) in figure 2 the reference plane shift is given by, 

Azp = -l p because A z° p is measured from a laminate’s original reference plane to the new reference 

plane, while l p is measured in the opposite direction, from the new reference plane to the laminate’s 
original reference plane. Considering first the standard laminate ABD matrix terms, 


A Group =Z A r = Z A J^ A Group = Z A " 


(53) 


Boroup = Z ®r = Z K + A °P ) : = Z K - ' l P A °P ) B Group = Z ( B J “ lp K ) 

k P P P 

Dcroup =Z D r=XK +2Az / ) K + V K ) 


(54) 


=ZK- 2 / p b ° +/ / a °) 


D 


Group 


=ZK- 


2/ p Bj+//A 


■J) ( 55 ) 


Similarly, the homogenized thermal terms are given by, 

^=Z(c) =ZK), 




(56) 
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Bg™,=Z(BS ew ) p =I[(B?) p+ ^(A?) p j 




n a 

^Group 


=e[K),-'*K ),5 

»§™p=ZK.) 

p p 

= ?[( DS l +2 *' ( BS l + *’ ( A? l^ S[(Bf ) # -2 /, (B?) p + V (A?) ( 

=l[K), - 2, pK) i> + '/K),,i! 

NLup=ZKew) =ZK) ^NLup=Z( N o)„ 

P p P 


MLup = Z ( M new ) = Z ( N 0 + ^ M 0 ) = Z ( N 0 - k ) “> M^up = Z ( M 0 “ h *0 ) 

P P P P P P P 1 

Finally, the homogenized thermo-electro-magnetic terms are given by, 

A aoup=Z( A o) ,j = E,M,ET,MT 


b L,p=Z|(bO -l p [Al)l,j=E,M,ET,MT 


•Cup'lN - 2 'p( b i) +'p 2 ( a 0 Ij=et,mt 


P ■" 


N Lup=Z( N o) ,j = E,M,ET,MT 


M Lup =Z( M 0 ->p N o) - i =E ’ M ’ ET - MT 


The homogenized constitutive equation for the group of laminates with respect to the reference plane 
defined by the origin of the z new coordinate in figure 2 is given by, 


( 57 ) 

1 (58) 

y 

(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 
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-J (66) 


4. Analysis of Blade Stiffened Panels 

The geometry of a blade stiffened panel is shown in figure 3. The face sheet can be an arbitrary 
composite laminate that is oriented as shown in the jc, y, z coordinate system such that the z-direction is 
the laminate through-thickness direction. This coordinate system is also used for the stiffened panel as a 
whole. The blade stiffener can also be an arbitrary laminate, but, as shown, it is oriented in the x B , y B , z B 
coordinate system with the blade through-thickness direction, z B , corresponding to the panel y-direction. 
Note that the x B -direction for the blade corresponds to the x-direction of the face sheet and panel as a 
whole. 

The key assumption made in incorporating the blade contribution to the overall behavior of the 
stiffened panel is that the contribution is largely decoupled from the face sheet contribution. That is, the 
response of the blade in its own y B -direction is assumed not to affect the panel and face sheet response 
due to the small contact region between the components. Considering first the panel mechanical behavior, 
the decoupling assumption is embodied by assigning the blade the following effective in-plane properties 
(in the local blade coordinate system), 

E?* 0, Ey =0, v*=0, G%= 0 (67) 

Thus, in terms of its effect on the panel response, the blade will have a contribution to the axial stiffness 
(x-direction), but no direct stiffness contribution and no Poisson contribution in the y B -direction (panel 

z-direction), and no shear contribution. E* is the effective (homogenized) axial stiffness of the blade. It 
can be calculated from the inverse of the blade laminate extensional stiffness matrix, A -1 , as, 



where t B is the thickness of the blade (see fig. 3) and Alf x is the 1 1 component of the inverse of the blade 

laminate extensional stiffness matrix. Based on the effective properties given in eq. (67), the uncoupled 
reduced stiffness matrix of the blade is, 

E* 0 0l 

Q u b = 0 0 0 (69) 

0 0 0j 
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Figure 3. — Geometry of a blade stiffened panel. The face sheet is shown in blue, 
while the blade stiffener is shown in pink. 

In the panel coordinate system (see fig. 3), the blade stiffness contribution is taken to be the volume- 
weighted sum of the blade uncoupled reduced stiffness matrix and the empty space between the blades 
along the panel y-direction. If the blade spacing (i.e., distance between adjacent blades) is denoted as S, 
this blade contribution is given by, 


qT 



e b 
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'1 
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0 
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0 

~SAll l 
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0 

0J 


(70) 


Employing eq. (2), the blade contributions to the ABD terms with respect to the blade’s midplane (see 
fig. 3) are given by, 


a» =or 


^ -h 3 h 5 - h 3 1 _ - cont 


» ° b =-\qT 


2 !) 
\2 


= QT(h 3-M 




= 0 


(71) 


(72) 


d°b =^Q“ nt 


h 5 - h 3 




*5-*3f 1 = ±Q r(A3 _ % )3 

J 12 


(73) 


Shifting these terms to the midplane of the face sheet, which serves as the reference plane for the panel, 

he -I - h'l 

using eqs. (9) to (1 1) with A z = — — yields, 
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Ag* =03- (%-*,) 


(74) 



Then, as discussed in section 3, the contributions of the blade and the face sheet to the panel ABD terms 
can be added to yield the effective panel ABD terms (since the blade and face sheet terms are now known 
with respect to the same reference plane). Thus, 


A panel ~ A FS + Q^”* (^3 “ % ) ( 77 ) 

Bpane^Bps+^Qr^S 2 -^) ( 78 ) 

Dpanel = »FS + \qT (*J "*?) ( 79 ) 


where A F s, B FS , and D FS are the face sheet laminate extensional, coupling, and bending stiffnesses, 
respectively. 

Considering the thermal response of the blade stiffened panel, the decoupling of the blade from the 
panel thermal response indicates that only the jc-direction thermal expansion of the blade will have an 
effect. Thus, the effective coefficients of thermal expansion (in the blade coordinate system, fig. 3) 
assigned to the blade for its thermal contribution are, 

a® ^ 0, a.y = 0, a^y = 0 (80) 


a® is the effective thermal expansion coefficient in the jc-direction. This term can be calculated by 
considering the inverted form of the general laminate constitutive eq. (47), 
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(81) 


Clearly, a thermal matrix can be identified as “CTE-like”, and, as identified by Collier (1993) (ref. 29), 
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B° 
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T>CC 

L^o 
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a lam 
coupling 
_ lam 


c- coupling 1 
°lam 

^lam J 


( 82 ) 


where a lam is a 3x1 vector relating strain to temperature change, af^ pling is a 3x1 vector relating 
curvature to temperature change, 5f^ pling is a 3x 1 vector relating strain to through-thickness temperature 
gradient, and 8i am is a 3x1 vector relating curvature to through-thickness temperature gradient. The 
x-direction effective thermal expansion coefficient for the laminate is then simply the first component of 
a lam , that is, 


a x = ( a lam)i 


(83) 


where (a lam ) ] is the 1,1 component of the 6x2 matrix formed by eq. (82), and the blade CTE vector is 


given by, 



a b = 


0 

°j 


(84) 


Employing eq. (1 5), the blade contributions to the panel thermal ABD terms with respect to the blade’s 
midplane (see fig. 3) are given by, 








-h 3 >2 

2 J 



T 1 

n 


= 0 


(85) 


(86) 




S* 


h ~^3 
v' 2 J 


( h. - h. s? ] 


- ^Qs nt a B (h “ % ) 3 


(87) 


Shifting these terms to the midplane of the face sheet, which serves as the reference plane for the panel, 
using eqs. (21), (23), and (27) with Az = — + — yields, 


( \Cont _ 

) = Qs° n a B (^3 - ^5 ) 



3 2 ~ Qs° nt a B (*3 “ *5 ) 



( 88 ) 

(89) 
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Qs Dt a B (^3 “ ^5 ) --Qs° nt a B (^3 - ^5 ) 


Then, as discussed in section 3, the contributions of the blade and the face sheet to the panel thermal ABD 
terms can be added to yield the effective panel thermal ABD terms (since the blade and face sheet terms 
are now known with respect to the same reference plane). Thus, 


A panel “ A FS + Qa”* a B (^3 - ^5 ) 

(91) 

® panel = ®FS + ^ Qa”* 3 “ ^5 ) 

(92) 

® panel = D FS + “Qa” 1 * a B (^3 _ ^5 ) 

(93) 


where A“ s , B“ s , and D“ s are the face sheet laminate thermal ABD matrices. 

A similar treatment applies to the thermo-electro-magnetic terms. The decoupling of the blade from 
the panel thermo-electro-magnetic response indicates that only the jc-direction electric, magnetic, thermo- 
electric, and thermo-magnetic expansion of the blade will have an effect. Thus, the effective piezoelectric, 
piezomagnetic, pyroelectric and pyromagnetic coefficients (in the blade coordinate system, fig. 3) 
assigned to the blade for its thermo-electro-magnetic contributions are, 

f —B -B -B 1 

e U e 2\ e 3\ 1 

0 0 0 

e fi = ° ° ° (94) 
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where e ; y , q t j , , and \|/j are the effective piezoelectric, piezomagnetic, pyroelectric, and 

pyromagnetic coefficients, respectively, that influence the jc-direction. In eq. (81) the following 
identifications can be made, 
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(98) 


(99) 


( 100 ) 


( 101 ) 


where e lam and q lam are 3x3 matrices relating strain to electric and magnetic field, respectively, 
e ton Plmg an ^ qiam Pl ' n8 are 3 X 3 matrices relating curvature to electric and magnetic field, respectively, 
ae lam and aq lam are 3x1 vectors relating strain to temperature changes (due to thermo-electric and 
thermo-magnetic coupling, respectively), aej ; °“ plmg and aqj^ pling are 3x1 matrices relating curvature to 
temperature change (due to thermo-electric and thermo-magnetic coupling, respectively), 5e f°^ p ling and 

gqcoupbng 

are 3x1 matrices relating strain to through-thickness temperature gradient (due to thermo- 
electric and thermo-magnetic coupling, respectively), and 8e lam and 8q lam are 3x1 matrices relating 
curvature to through-thickness temperature gradient (due to thermo-electric and thermo-magnetic 
coupling, respectively). 

As was the case for the blade thermal expansion coefficients, the thermo-electro-magnetic 
coefficients present in eqs. (94) to (97) can be extracted from eqs. (98) to (101). The ejj and qfj terms 

are simply the first rows of the e lam matrix and the q lam matrix, respectively, while , and \j/f are 
simply the first term in the ae lam and aq lam vectors. 

Employing eqs. (36), (38), (45), and (46), the blade contributions to the panel electric, magnetic, 
thermo-electric, and thermo-magnetic ABD terms with respect to the blade’s midplane (see fig. 3) are 
given by, 
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Shifting these terms to the midplane of the face sheet, which serves as the reference plane for the panel, 

he “I" h'l 

using eqs. (48, 49) with A z = — — yields, 


( A f ) C ° nt = Qa° nt e g (h 3 - 1 %) 
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(^B ) _ Qs° nt <Js (^3 - ^5 ) (1 14) 

(b % ) COnt =^^QT q b {h 3 - h 5 ) = ^QT q b [h 3 -hi) (115) 

( Af 7 ) C ° nt = Q^ ont l B (h 3 - h 5 ) (116) 

( B f ) C ° nt =^QT Ib (h 3 -h 5 ) = ^QT Ib ~hi) (117) 

(Df 7 ) C ° nt =^QT Ib^- h 5 f + j Qg 0nt l B {h 3 - h 5 )=^QT «b (^3 - h] ) (118) 

(Arf° nt =Q Tv B {h 3 ~h 5 ) (119) 

(Bf J 00 " = ^Q r Vb (h 3 ~h 5 ) = \QT Vb (hi -hi) (120) 


(i2i) 

Then, as discussed in section 3, the contributions of the blade and the face sheet to the panel electric, 
magnetic, thermo-electric, and thermo-magnetic ABD terms can be added to yield the corresponding 
effective panel ABD terms (since the blade and face sheet terms are now known with respect to the same 
reference plane). Thus, 


panel = A FS + Qb”* e B {h 3 ~ h 5 ) 

(122) 

Bpanel = ®FS + ^QB° nt e B (^3 “^5 ) 

(123) 

A panel = ^-FS + Qb^ q« (A 3 ~ h 5 ) 

(124) 

Bpanel = B fs + ~ Qb”* qg (h% " ^5 ) 

(125) 

bpanel = ^FS + Qb”* ^B (h 3 ~ h$ ) 

(126) 

Bpanel = B fs + ~ Qa”* £s (^3 ~ h$ ) 

(127) 
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® panel “ D FS + J Qb”* (*3 “ ^5 ) 

(128) 

a£Li = A^ +QS* ¥*(*-*) 

(129) 

® panel = B Fsf + ^Qs° nt Vfi (^3 “ ^5 ) 

(130) 

D panel = D FsT + “Qs”” 1 (^3 “ ^5 ) 

(131) 


where Ap S , Bp S , and Dp S are the face sheet laminate electric, magnetic, thermo-electric, and thermo- 
magnetic ABD matrices. 

It should be noted that if the blade stiffened panel includes a bottom face sheet, its contribution can be 
included by simply shifting its mechanical, thermal, electric, magnetic, thermo-electric, and thermo- 
magnetic ABD terms to the face sheet midplane according to eqs. (9) to (1 1), (21), (23), (27), (48), and 
(49) and adding these shifted terms to the corresponding panel ABD terms. 


5. Analysis of Panels with Flanged Stiffeners 

Considering a flanged stiffened panel, as shown in figure 4, the outlined region is referred to as 
“Segment 8”. Comparing figure 4 to figure 3, it is clear that the only difference between the flanged 
stiffened panel and the blade stiffened panel is the presence of the flange. Thus, the development of 
effective ABD terms for the flanged stiffened panel necessitates only the additional analysis of Segment 
8. As shown in figure 5, Segment 8 is composed of two other segments, Segments 1 and 6. 

In the ^-direction, which is the direction of the stiffener (and thus the strong direction), iso-strain and 
iso-curvature conditions are employed. It is further assumed that the homogenized Segment 8 jc-direction 
force and moment resultants are the sum of the volume- weighted jc-direction force and moment resultants 
of the constituent segments (Segment 1 and 6). This is analogous to the “mechanics of materials” 
micromechanics approach presented by Herakovich (ref. 28) that homogenizes continuous fiber 
composites to obtain effective properties. The conditions employed for the stiffener direction are, 


(4M4) 6 «(4) b 

( K «) 1= ( K xr) 6 =(Kxc) 8 

( *« f = j~ [{FWNT) (N„ f + ( WNT ) (JVjpp ) 6 1 
(^x*) 8 =^-[(FB7VT)(M JCC ) 1 + (WNT)(M xx f 


(132) 

(133) 

(134) 

(135) 


where Sx is the stiffener spacing, FWNT is the distance between flanges, and WNT is the flange width. 
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Figure 4. — A flanged stiffened panel with the segment 8 region outlined. 
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Figure 5. — The Segment 8 region divided into Segments 1 and 6. 


In the y-direction and for the in plane shear effect (xy-components), the converse assumptions are 
employed. That is, iso-force and moment resultant conditions are imposed, along with volume-weighted 
summing for the strain and curvature components. These conditions are, 
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The constitutive equation for each segment is given by, 
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where B T is the transpose of the B matrix, which is necessary because the B matrix itself is typically 
asymmetric for Segment 8 (although the ABD matrix remains symmetric). The equations are rearranged 
such that all components that fall under the iso-assumptions appear on the left hand side, 
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(145) 


or, in simplified notation, 
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(146) 


where the hatted matrices represent ABD matrices altered via the rearrangement required in forming eq. 
(145) from eq. (144). Solving eq. (146) for R B yields, 


NASA/CR— 2009-21 5269 


26 



( 147 ) 



Now, employing the above simplified notation, all of the iso-assumptions can be expressed as, 

(Ra) 1= (Ra) 6 =(Ra) 8 (148) 

while all of the volume-weighted summation assumptions can be expressed as, 

(R b ) 8 =-^-[(Flf7V7’)(R B ) 1 +(lf7Vr)(R B ) 6 J (149) 


The R b vector for each segment in eq. (149) can be replaced with the expression from eq. (147), 



Each effect present in eq. (150) can be isolated by considering the case when only it is active. That is, if 
we assume that the panel has only mechanical loading and no thermal, electric, or magnetic loading, 
eq. (150) simplifies to, 

S* [C 8 (R A ) = (FWNT) [C 1 ] _1 (R A ) + (WNT) [c 6 ]"' (R A ) (151) 
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A Q 

Then, inverting eq. (151) provides C , and rearranging the terms (see below) provides the homogenized 
ABD matrix for Segment 8. 

For the thermal, electric, magnetic, thermo-electric, and thermo-magnetic terms appearing in 
eq. (150), we have, 

Thermal 

^C 8 ]- 1 [C« 8 ] ^^|j= (F’WW7’)[c 1 * C al ^J j + (WNT)[c^[c^] A ^|j (153) 

or 

c“ 8 =^- C 8 ^(FlF?V7 , )[c 1 -1 C al ] + (lF?V7 , )[c 6 C a6 | (154) 

Electric 

i r*i i w -i ^ 

^[c 8 ] [c £8 ] E 2 C- E 2 + (1F7V7 , )[c 6 ] C- E 2 (155) 

_e 3 j e 3 j e 3 j 

or 

C £8 ]=-L c 8 |(F0W7 , )[c 1 J 1 [c £1 ] + (fF?V7 , )[c 6 p 1 [c £6 ]J^ (156) 

Magnetic 

r*.i r^ii r^ii 

^[c 8 ] [C-] H 2 = (FW2V7’)[d 1 ] [c M1 ] H 2 +(lfwr)[c 6 ] [c^ 6 ] H 2 (157) 
h 3 j h 3 j h 3 j 

or 

C M8 ] = y- C 8 ] [(FWNT^J 1 [c“'] + (WNT)[c*J l [c- 6 ]| (158) 

Thermoelectric 

+ (JWr) [c* ]“' [c CT6 ] _^j (159) 

or 

[c ET « =^[c 8 ]5(Fff7V7’)[c 1 ] _1 [c £n ] + (lf7Vr)[c 6 p 1 [c £7 ’ 6 ]| (160) 
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By rearranging each of eqs. (152), (154), (156), (158), (160), and (162), the Segment 8 thermal, 
electric, magnetic, thermoelectric, and thermomagnetic ABD matrices can be determined. However, for 
this to be accomplished, all of the hatted matrices for each segment must be determined. This is done by 
separating the segment constitutive equation as, 


N 


xy 


M 


yy 


M 


xy 


’22 

^23 


a 2 

5 2 “ 

B? 


’23 

4 33 


B 


’22 


■®23 ,1 


B?, 2 B 


33 


B22 B'xi D’ 


'32 


22 


'23 


\_ B 23 B 33 D 32 D 33\ 




yy 


+ 


M A 


*12 

4 13 


By\ ,1 


B 


31 


B\ 2 B) 


12 


r e ° i 

^xx 

_^xx J 


5 2 “] 

5“ 

d 3 «u 


A7 7 1 
-ACT] ' 


,£ 

A 2\ 


A E 

A 22 


a E “I 
^23' 1 


Bi ft B> 


l 32 

£ 


J 21 

:>£ 


*31 A 


'22 

£ 




'32 


,££ 

4 3 

d£^ 


D 


,££ 

2 

,££ 

'3 

,££ 


1 


ET 


Df\\ 


AT 1 
L-AG 1 !' 


aMT 


23 

^IJ 


A 3 J 




L^y 


jM 

A n 

a m 

a 3\ 

r>M 


BSt B 


'2\ 

r>M 


aMT 


B 


>MT 


) MT 


) MT 


d? t 

dT\\ 


AT 1 
L-AG 1 ] 


aM 

A 22 

a m 

A 32 
M 


*31 * 


'32 


a M 1 
^23 I 1 


22 

M 


^33 
B 23 I 

*33lJ 


'^ll 
^2 
1*3 J 


(163) 


or, employing simplified notation, 


v, = S, lUj + S 12 u 2 - S“AT - Sf E - sf H - Sf Y AT - Sf AT 


,ET, 


?MT , 


(164) 


and 
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or, employing simplified notation, 

v 2 = S 21 Uj + S 22 u 2 - S 2 AT - Sf E - Sf H - Sf 7 AT - S^AT 
Equations (164) and (166) can be rearranged and combined as, 
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(167) 


where now all components subjected to the iso-conditions appear on the left hand side. Solving eq. (166) 
for u 2 yields, 


u 2 = [S 22 ] 1 {v 2 - S 21 u, + S 2 AT + sf E + sf H + Sf r AT + S^AT} 
and substituting for u 2 in eq. (164) using eq. (168) yields, 

Vi=S llUl + S 12 [S 22 r 1 {v 2 - S 2lUl + S 2 AT + sf E + sf H + Sf 7 AT + Sf T AT] 


(168) 


- S“AT - Sf E - Sf H - Sf AT - Sf" AT 


?ET A 


i 


(169) 


or, with terms grouped, 

v i =|sii -s 12 [s 22 ] s 21 Ju 1 + s 12 [s 22 ] v 2 -{sf -s 12 [s 22 ] s 2 Jat 

-{sf - S 12 [S 22 r 1 sf } E - {sf - S 12 [S 22 r 1 sf } H - {sf 7 - S 12 [S 22 f 1 Sf } AT 
-{sf 7 -S 12 [S 22 ] _1 Sf 7 } AT 


(170) 


Comparing eqs. (170) and (168) with eq. (167) allows the identification of the following hatted terms, 

Mechanical 

1 = Si 1 - S 12 [S22 F 1 s 21 (171) 
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Thermal 


Electric 


Magnetic 


Thermoelectric 


Thermomagnetic 
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The equivalence between the S* matrix terms in eqs. (171) to (184) and the C* matrix terms from 
eqs. (152), (154), (156), (158), (160), and (162) is, 
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(190) 

where the parenthetical indices refer to the indices of the components within each S t ' 

and S* matrix. 

Using eqs. (152), (154), (156), (158), (160), and (162), we can now determine the hatted quantities for 
Segment 8 from the hatted quantities of Segments 1 and 6 given by eqs. (171) to (190). In order to then 
determine the non-hatted quantities for Segment 8, the reverse of the procedure embodied by eqs. (164) to 
(190) are employed, resulting in equations of identical form of eqs. (171) to (184), but with the rolls of 
hatted and unhatted quantities reversed. That is, 
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( 200 ) 

( 201 ) 

( 202 ) 

(203) 

(204) 


Finally, the form of the equivalence between the S matrix terms and the C matrix terms is identical to the 
form of the equivalence between the S matrix terms and the C matrix terms given in eqs. (185) to (190). 


6. The Thermal Analogy 

It is clear from the preceding development that an analogy exists between the thermal treatment of the 
material/laminate and each of the piezoelectric, piezomagnetic, thermoelectric, and thermomagnetic 
treatments of the material laminate. That is, the piezoelectric, piezomagnetic, thermoelectric, and 
thermomagnetic effects are each incorporated into the overall theory in a way that is analogous to the 
treatment of thermal effects. Therefore, it is possible to utilize the thermo-electro-magnetic terms within 
the theory to mimic thermal behavior, and conversely, it is possible to utilize the thermal terms within the 
theory to mimic the piezoelectric, piezomagnetic, thermoelectric, and thermomagnetic effects. This 
thermal analogy can be used to verily the thermo-electro-magnetic implementation vs. an established 
thermal implementation or to perform a piezoelectric, piezomagnetic, thermoelectric, and thermomagnetic 
analysis using an established thermal implementation. Cote et al. (ref. 30) employed the latter approach to 
simulate the dynamic response of composite beam with an embedded piezoelectric actuator. 

6.1 Material Level Thermo-Electro-Magnetic Terms Mimicking Thermal Response 

The general thermo-piezo-electro-magnetic constitutive equation for a material is given by, 
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This equation can be written in matrix form as, 

o = C [e - a AT] - e [E + 5 AT] - q [H + y AT] 
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Thus, in order to mimic the material level thermal effects using the piezoelectric terms, we seek a ' 
piezoelectric coefficient matrix such that a AT = C -, e 
Assuming orthotropic behavior, 
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Now, by setting E ] = AT, E 2 =0, E 3 = 0 , we have, 
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( 209 ) 


( 210 ) 


( 211 ) 


Thus, by utilizing the fake piezoelectric coefficients indicated by eq. (211), along with E x = AT, E 2 = 0, 

E 3 = 0 (and then also utilizing fake thermal expansion coefficients, a Fake = 0 ), it is possible to determine 
the thermal behavior of a material through the piezoelectric terms. This can be used to verify the 
piezoelectric implementation. 

Similarly, to verify the piezomagnetic effects, we set H x = AT, H 2 =0, H 3 = 0 , and 
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For the thermo-electric effects, what we have a AT = C 1 e Fake £ Fake AT . Again, assuming 
orthotropic material behavior, 
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Setting =1, ^2 ake =0, Cs ake =0, we have 
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Similarly, for the thermo-magnetic terms, we set \|/j =1, \|/ 2 = 0, \|/ 3 = 0 and obtain, 
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6.2 Material Level Thermal Terms Mimicking Thermo-Electro-Magnetic Response 

In order to mimic the material piezoelectric behavior using the thermal terms, we again begin with 
eqs. (205) to (207). Setting a Fake AT Fake = C -1 eE and retaining anisotropic behavior, 
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Setting A7 Fake = E ] , E 2 =0, E 3 = 0 , we have, 
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Generally, for any single applied electric field component, E t , set Ar Fake = E i and Ej = 0 (y ± i) along 
with, 
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Similarly, to mimic the piezomagnetic effects using the thermal terms, set A7’ Fake = H i and H • = 0 
( j i ) along with, 
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In order to mimic the thermo-electric effects, we set a Fake A7’ Fake = C ] e^AT . Retaining anisotropic 
material behavior and setting A7’ Fake = AT , 
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In order to mimic the thermo-magnetic effects, we set Ar^ 6 = AT , along with, 
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It is thus possible to mimic the piezoelectric or piezomagnetic behavior of a material in response to a 
single electric or magnetic field component by using the fake coefficients of thermal expansion given by 
eqs. (221) and (222). The thermo-electric and thermo-magnetic behavior can be mimicked by using the 
fake coefficients of thermal expansion given by eqs. (223) and (224). 

6.3 Panel Level Thermal Terms Mimicking Thermo-Electro-Magnetic Response 

A simple thermal analogy is also in effect on the panel and laminate level. The laminate or stiffened 
panel constitutive equation is given by, 
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Equating fake thermal terms with the piezoelectric terms gives, 
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or, 
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Applying a single electric field component, E h we have, 
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Setting A7 Fake = E t and AG Fake = 0 gives, 
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Similarly, to determine fake thermal terms that mimic the panel/laminate piezomagnetic response, we set 
AyFake _ an( j AG Fake = 0 along with, 


4 

-jFake 

a m~\ 

A u 

4 


a m 

A 2 i 

4 


a m 

Ai 



r m 

b m 

4 


t>M 

B 2 i 

4 

J 

n M i 
L°3i J 


( 231 ) 


NASA/CR— 2009-21 5269 


41 



These results are useful because they enable the analysis of piezoelectric and piezomagnetic shells 
within NASTRAN using the software’s thermal analysis capabilities. NASTRAN accepts thermal ABD 
matrices for the shell materials. Thus, by providing NASTRAN with the appropriate fake thermal terms, 
the software will solve a thermal problem that is analogous to a desired piezoelectric or piezomagnetic 
problem. An even simpler analogy exists between the panel level thermal terms and the thermo-electric or 
thermo-magnetic terms. Simply setting, 
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will enable the thermal capabilities to mimic the thermo-electric or thermo-magnetic response. However, 
these analogies are less useful because their use assumes that the panel or laminate only reacts to the 
applied thermal loading through the thermo-electric or thermo-magnetic terms, and not through the 
standard thermal expansion terms. That is, these analogies would only be valid for the case in which the 
panel/laminate coefficients of thermal expansion are zero. 

7. Results and Discussion 

To verify the thermo-electro-magneto-elastic implementation within HyperSizer, we consider a 
bonded facesheet-flange combination, which is identical to the HyperSizer section 8 shown in figures 4 
and 5. The geometry of this problem is shown in figure 6. Within HyperSizer, this case corresponds to a 
T-stiffened panel (see fig. 4) with an infinitesimal web. The facesheet is composed of PZT-7A zirconium 
lead titanate piezoelectric material with a through-thickness poling direction, while the flange is 
composed of aluminum. The material properties employed for these materials are given in tables 1 and 2. 
As indicated in figure 6, a voltage difference of 1 x 10 5 V is applied through the thickness of the face sheet. 
This corresponds to an electric field of (1*10 5 V)/(0.00229 m) = 43.67 MV/m. This problem has also 
been analyzed using the ABAQUS finite element analysis package, employing the finite element mesh 
shown in figure 7, consisting of a total of 21,200 elements. The facesheet is composed of CPE4E 
piezoelectric plane strain continuum elements, while the flange and adhesive are composed of CPE4R 
reduced integration plane strain continuum elements. 
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Figure 6. — ABAQUS solution domain for a facesheet-adhesive-flange combination. 
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TABLE 1— ELASTIC CONSTANTS OF THE MATERIALS USED IN THE ANALYSES (REFS. 3 AND 4) 



E„ 

(GPa) 

E 22 

(GPa) 

E 33 

(GPa) 

V 12 

V 13 

V 23 

G 12 

(GPa) 

G 13 

(GPa) 

g 2 3 

(GPa) 

Aluminum 

PZT-7A 

68.95 

94.97 

68.95 

81.90 

68.95 

94.97 

0.30 

0.384 

0.30 

0.323 

0.30 

0.331 

26.52 

25.40 

26.52 

25.40 

26.52 

35.90 


TABLE 2,— PIEZOELECTRIC PROPERTIES (REF. 4) OF THE MATERIALS USED IN THE ANALYSES 



e 222 

(C/m 2 ) 

®233 

(C/m 2 ) 

e 211 

(C/m 2 ) 

e 323 

(C/m 2 ) 

e 112 

(C/m 2 ) 

k?? 

(10 -9 C/V m) 

^33 

(1CT 9 C/Vm) 

kn 

(KT 9 C/V m) 

PZT-7A 

12.25 

-2.1 

-2.1 

9.2 

9.2 

2.07 

4.07 

4.07 


The normal stresses in the plane defined by the facesheet (x-y plane) arising due to the applied voltage 
along the cut shown in figure 7 are plotted in figures 8 and 9. The original ABAQUS solution (labeled 
“ABAQUS Fully Coupled”) agrees reasonably well with the HyperSizer solution, but some deviation is 
evident. This deviation is caused by the fact that HyperSizer’s formulation is based on the ability to apply 
constant panel (or laminate) level electric field components. In contrast, the ABAQUS continuum 
solution involves the application of electric potential at boundaries along with solution of a boundary 
value problem for the electric potential throughout the model. As will be shown, the ABAQUS solution 
does not result in constant electric field components. 
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Figure 8. — Stress component (a x ) along the cut through the flange adhesive 
and facesheet defined in figure 7 as predicted by HyperSizer and 
ABAQUS both with fully coupled and decoupled electric field components. 
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Figure 9. — Stress component (a y ) along the cut through the flange 
adhesive and facesheet defined in figure 7 as predicted by HyperSizer 
and ABAQUS both with fully coupled and decoupled electric field 
components. 
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In order to provide the ABAQUS solution with constant electric field components, which mimics the 
HyperSizer formulation, it is possible to employ artificially large dielectric constants (1%) for the PZT-7A 
material. This, in effect, decouples the electric field components (E t ) in the piezoelectric material 
constitutive equation, 

°lll C U Q 2 Q 3 0 ^ 0 e ll 0 0 1 Bnl 

°22 Q2 ^22 ^23 ® 0 0 e 12 0 0 8 2 2 

O33 C 13 C 23 C 33 0 0 0 e 13 0 0 e 33 

a 23 0 0 0 C44 0 0 0 0 0 y 2 3 

a 13 = 0 0 0 0 C55 0 0 0 e 35 y 13 (234) 

a i2 0 0 0 0 0 C 66 0 e 2 6 0 Y12 

Dj I 6]2 C13 0 0 0 Kj j 0 0 £) 

D 2 0 0 0 0 0 e 2 6 0 k 22 0 E 2 

D 3 j [ 0 0 0 0 e 35 0 0 0 k 33 J |_ £3 J 

where D t are the electric displacement components and an jci poling direction has been assumed. By 
making the dielectric constants large, the stress, strain, and electric displacement components will depend 
on the electric field components, but the electric field components will not depend on the stress, strain, or 
electric displacement components. This condition is analogous to specifying constant electric field 
components, as is done in HyperSizer, such that they cannot vary due to the other field components. 

To accomplish the decoupling of the electric field components described above, the dielectric 
constants given in table 2 were increased by a factor of 10 5 . Results for this case, labeled “ABAQUS 
Decoupled”, in figures 8 and 9 now agree extremely well with the HyperSizer solution. Figures 10 to 13 
provide further comparisons between the ABAQUS results with and without fully electric field coupling. 
Figure 10 indicates a slight difference in the von Mises stress field between the two cases. Likewise 
figure 1 1 shows only a small difference between the electric potential solution between the fully coupled 
and decoupled cases. However, the spatial derivatives of the electric potential, which are the electric field 
components in the two directions, do show significant differences. In figure 12(a), in the region of the 
facesheet directly beneath the flange, a through-thickness electric field gradient has arisen due to stress, 
strain, and electric displacement components that arise in the region due to the presence of the flange. In 
figure 12(b) on the other hand, the decoupling has eliminated this variation in the E z electric field 
component, and a constant E z value of 43.67 MV/m results in the facesheet, identical to the value applied 
in HyperSizer. The variation in E z evident in figure 12(a) is approximately ±8 percent with respect to the 
constant decoupled value, which is not excessively large, but clearly is large enough to have an effect. In 
figure 13(a), a gradient in the E x electric field component has arisen near the free edge of the adhesive 
bond within the facesheet. This is due to y xz shear strain that arises in this region which, thanks to a non- 
zero e 2 6 (see eq. (234)), gives rise to an electric displacement and field. In figure 13(b), the decoupling of 
the electric field components has eliminated this gradient, and a constant E x field results with a value of 
zero. Again, this is identical to the condition imposed in the HyperSizer solution. 
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Figure 10. — Von Mises stress field predicted by ABAQUS with 
(a) fully coupled and (b) decoupled electric field components. 
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Figure 1 1 . — Electric potential field predicted by ABAQUS with 
(a) fully coupled and (b) decoupled electric field components. 
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Figure 12. — Electric field component E z predicted by ABAQUS with 
(a) fully coupled and (b) decoupled electric field components. 
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Figure 13. — Electric field component E x predicted by ABAQUS with 
(a) fully coupled and (b) decoupled electric field components. 
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8. Conclusion 


The methods employed within HyperSizer to analyze composite stiffened panels with thermo-electro- 
magneto-elastic coupling have been presented. Starting on the level of the laminate, classical lamination 
theory is employed, and new electric, magnetic, thermo-electric, and thermo-magnetic ABD terms have 
been identified. Homogenization techniques have been presented for blade and flanged stiffeners that 
result in stiffened panel level constitutive equations that are analogous to thermo-electro-magnetic 
laminate constitutive equations. These allow the stiffened panel to be included in a higher scale structural 
model via methods that accept laminate constitutive equation terms such as MSC/NASTRAN. While 
MSC/NASTRAN does not presently accept the newly developed thermo-electro-magnetic terms, the 
ability to calculate these terms allows the easy quantification and analysis of a new class of smart 
stiffened panels. 

The analogy between standard lamination theory thermal expansion effects and the newly developed 
thermo-electro-magnetic effect was also discussed. This thermal analogy can be used to model the 
response of smart structure through a method’s existing thermal capabilities, or to model a thermal 
response of a structure through a method’s electro-magnetic material capabilities. The former is useful for 
generating results for comparison with a known method, while the latter is useful for verifying a thermo- 
electro-magnetic method by generating known thermal results. In both cases, the coefficients needed to 
take advantage of the thermal analogy have been identified. 

Verification results have been presented that compare HyperSizer piezoelectric results with ABAQUS 
piezoelectric finite results for a facesheet-adhesive-flange combination. Because the ABAQUS 
capabilities are limited to piezoelectric materials, the new piezomagnetic, thermo-electric, and thermo- 
magnetic capabilities of HyperSizer were not compared. The results indicated generally good agreement 
between HyperSizer and ABAQUS. However, because ABAQUS solves a piezoelectric boundary value 
problem subject to prescribed electric potential, the electric field components can vary spatially within the 
model. The HyperSizer implementation, on the other hand, is based on classical lamination theory with 
prescribed spatially constant electric field components. This limitation leads to some discrepancy in the 
results. By decoupling the electric field components within the ABAQUS solution by significantly 
increasing the piezoelectric material dielectric constants, the electric field components are forced to 
remain constant. ABAQUS results in this decoupled condition match the HyperSizer results nearly 
exactly, serving to verify the HyperSizer piezoelectric implementation. 
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